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ABSTRACT 

Winds from accretion disks have been proposed as the driving source for precessing 
jets and extreme bipolar morphologies in Planetary Nebulae (PNe) and proto-PNe 
(pPNe). Here we apply MHD disk wind models to PNe and pPNe by estimating 
separately the asymptotic MHD wind velocities and mass loss rates. We discuss 
conditions which may occur in PNe and pPNe accretion disks that form via binary 
interactions. We show that the resulting winds can recover the observed momentum 
and energy input rates for PNe and pPNe. High accretion rates (M a ~ 10 M yr _1 ) 
may be required in the latter case. We find that the observed total energy and 
momentum in pPNe can be recovered with disk wind models using existing disk 
formation scenarios. When combined with existing scenarios for accretion disk 
formation from disrupted stellar companions, our models may provide an explanation 
for the existence of high speed polar knots (FLIERS) observed in some PNe. 



Subject headings: Planetary Nebulae — magnetic fields — magnetohydrodynamics: 
MHD — shock waves 



1. INTRODUCTION 



1.1. Introduction 



Planetary Nebulae (PNe) and proto-Planetary Nebulae (pPNe) are believed to be the 
penultimate evolutionary stages of low and intermediate mass stars (M* < 8 M Q ). PNe and 
pPNe appear on the sky as expanding plasma clouds surrounding a luminous central star. As 
ground based telescopes increased their resolution, elliptical and bipolar PNe were revealed (Balick 
(1987)). The bipolar nebulae may be further subdivided into "butterfly," shapes in which the 
"waist" is pinched at the central star, and "bilobed" PNe in which a pair of larger outer lobes 
connects to a central and generally smaller round or elliptical nebula (for a review see Balick & 
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Frank (2002) and references therein). More recently, deeper and higher resolution studies with 
the HST have shown evidence for narrow collimated features that appear to be better described 
as jets than bipolar lobes. The nature of these jets and other highly collimated bipolar outflows in 
both PNe and pPNe remains a subject of considerable debate (Sahai & Trauger (1998)). 

Considerable progress has been made in understanding the hydrodynamic shaping of elliptical 
and bilobed PNe (for reviews Frank (1999),Balick & Frank (2002)), however the origin of extreme 
butterfly nebulae as well as jets in PNe continues to pose a number of problems for theorists. The 
outstanding issues include: (1) lack of a large scale torus to focus the outflows, (2) the very large 
velocities in some jets, (3) point-symmetry of outflows. 

Beyond issues of morphology, there exists a formidable problem for pPNe concerning the 
total momentum and energy in the outflows. A number of observational studies have shown 
that radiatively accelerated winds in pPNe cannot account for the high momentum and energy 
implied by CO profiles. This problem was identified first by Knapp (1986), and most recently and 
comprehensively investigated by Bujarrabal, Alcolea and collaborators (Bujarrabal et a/.(2001)). 
The latter find that 21 of 23 CO-emitting pPNe objects show outflows whose scalar momentum 
(n = MV) are more than 10 3 times larger than that in the stellar radiation. Thus n >> (L*/c)At, 
where here L* is the stellar luminosity emitted during the pPNe outflow expansion lifetime At. In 
light of these results both the launching and collimation of winds in pPNe becomes problematic. 

The dominant hydrodynamic theory for shaping PNe had been the Generalized Interacting 
Stellar Winds (GISW) model (Kwok et a/.(1978), Balick (1987), Icke (1988)) in which a star 
and its wind evolve from the AGB to a white dwarf. A slow, dense (wind expelled during the 
AGB is followed by a fast, tenuous wind driven off the contracting proto-white dwarf during 
the PNe phase. Numerical models have shown this paradigm can embrace a wide variety of 
nebular morphologies including highly collimated jets (Icke et al. (1992), Mellema k, Frank(1997), 
Borkowski et aZ.(1997)) when the slow wind takes on aspherical density distributions. While the 
GISW model can produce narrow jets it usually requires a large-scale "fat" torus. It is difficult 
to imagine that a large-scale out-flowing gaseous torus can provide a stiff precessing nozzle for 
production of point symmetric flows (Recent results of Icke 2003 however indicate that pure 
hydrodynamics may lead to point symmetric shapes in some cases). 

In addition, it is now recognized that fast (> 100 km s _1 ) bipolar outflows can occur in the 
pPNe or even the Post- AGB stage. Objects like CRL 2688 and OH231. 8+4.2 raise the question of 
how high-velocity collimated flows occur when the star is still in a cool giant or even supergiant 
stage (CRL 2688 has an F Supergiant spectral type). Finally the GISW model assumes a radiation 
driven wind. As discussed above, the results of Bujarrabal et al. (2001) make radiation driving 
untenable as a source for many pPNe flows. 

Models invoking a toroidal magnetic field embedded in a normal radiation driven stellar 
wind have shown considerable promise (Chevalier & Luo (1994), Rozyczka, &: Franco 1996, 
Garcia-Segura et aZ.(1999)), Garcia-Segura (1997)). Recent results, however also imply that jets 
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may form at smaller distances and that such models may not begin with appropriate initial 
conditions (Gardiner & Frank (2001)). In addition, by their very nature, MWB models can not 
account for the momentum excesses in pPNe since they also require radiation driven winds. The 
fields are simply too weak to power the observed outflows. 

Thus there remains considerable uncertainty about the processes which produce collimated 
jets in pPNe and PNe. Other systems which produce jets such as YSOs, AGN and micro-quasars 
have been modelled via a combination of magnetic and centrifugal forces from accretion disks 
(Blandford & Payne 1982, Shu et a/. (1994), Konigl & Pudritz (2000)) 

The success of these Magneto-centrifugal Launching (MCL) models is such that it is 
worthwhile considering if such a scenario can be applied to PNe and pPNe. Indeed, Morris (1987) 
and Soker & Livio (1994) mapped out scenarios in which accretion disks form around binary PNe 
progenitors. Each study equated the existence of disks with the existence of jets. The details of 
the jet launching and collimation mechanism where not, however, specified. Recent work by Soker 
& Rapport 2001, Soker 2001 and Soker &, Livio (2001) have relied heavily on collimated winds 
from disks but these works also do not specify how such winds are launched or collimated. Thus 
application of MCL disk wind models to PNe and pPNe would close an important gap in building 
a new paradigm for these systems. In particular MCL models may offer a means of resolving issues 
associated with both PNe jet precession (the underlying disk precesses via instabilities: Livio & 
Pringle (1996), Quillen (2001)) and those associated with pPNe momentum excesses (no need for 
radiation pressure launching winds). 

Recent studies by Blackman, Welch & Frank (2001) and Blackman et al. 2001 have explored 
MCL paradigms for both the star and disk in pPNe and PNe systems. They computed an upper 
limit on the magnetic luminosity available to power an outflow, assuming that a dynamo is the 
source of the large scale magnetic field. In this paper we we provide further calculations along 
these lines deriving scaling relations from the equations for MCL and separately estimate the mass 
outflow rate and the asymptotic outflow velocity. We then compare to these results to PNe and 
recent observations of PNe. When the dynamo is invoked to produce the field, the mechanical 
wind luminosity and thus outflow rate are naturally linked to the accretion rate. 

In section 2 we derive the outflow speed and the mass loss rate by combining results from 
magnetically driven wind theory and dynamos. In section 3 we discuss models of pPNe and PNe 
disks. In section 4 we apply the results of sections 2 and 3 to PNe and pPNe and show that it is 
easy to solve the afore mentioned the pPNe momentum and energy excess problems. In section 5 
we conclude and discuss open questions. 



-4 - 



2. Mass Outflow Rate and Asymptotic Wind Speed from MCL Theory and 

Dynamo Theory 

2.1. Magnetic Luminosity 

The basic physics of magneto-centrifugal launching of winds and jets is well studied when 
a magnetic field distribution is imposed on the disk (Blandford & Payne 1982, Sakurai (1985), 
Pelletier & Pudritz (1992) (hereafter PP92) Shu et al. (1994) Ostriker (1997), Konigl & Pudritz 
(2000) (hereafter KP00)). These models presume an initial field of a given strength and geometry, 
but the tendency for large scale fields to diffuse (Lubow & Pringle 1994; Blackman 2003) suggest 
that the field must be generated in situ by a dynamo (e.g. Blackman et al 1999). In this section 
we briefly summarize how to combine the basics of Poynting flux driven outflows with asymptotic 
wind solutions and mean field dynamo theory to estimate the asymptotic wind speed and the 
outflow accretion rate. 

Magneto-centrifugal launching is a means of converting gravitational binding energy in an 
accreting source into kinetic energy of an outflowing wind. The magnetic fields act as a drive belt 
to extract angular momentum from the anchoring rotator and launch the wind. The magnetic 
luminosity, or equivalently, the maximum magnetic power available for a wind can be obtained 
from the integrated Poynting flux (Blackman, Frank &: Welch (2001)). The field lines rotate 
nearly rigidly at angular speeds associated with the anchoring foot point ^o( r ) U P to the Alfven 
radius r^(r). After this point the angular speed falls off with 1/r 2 (conserving specific angular 
momentum) and the field falls off as 1/r. The field is primarily poloidal out to the Alfven radius 
associated with each field line where the poloidal and toroidal components are comparable. If the 
poloidal field falls off as 1/r 2 out to the Alfven radius and is not too far from the disk surface, the 
Poynting flux can be approximated by the contribution from the Alfven radius associated with 
field lines anchored at the inner radius of the disk. That is, we have 

1 f ["^AiTi) 

L w = -M^L ~ L mag = J (E x B) • dS A ~ J n(r)rB p B^rdr ~ B 2 A n r\, (1) 

where ro is the disk inner radius and where Ba = B^ ~ B p at the Alfven surface (i.e. the toroidal 
and poloidal field components are nearly equal). Since the dominant contribution to the magnetic 
luminosity comes from the field lines anchored at the innermost radius r,, in what follows all 
quantities labeled with subscript refer to those values evaluated at r = T{, the inner most disk 
radius. This also implies that ro = n. 

2.2. The Bernouilli Constant 

To estimate the mass outflow rate and outflow speed separately a bit more work is 
required. The MCL problem requires the construction of solutions for a steady, ideal, isothermal 
magnetohydrodynamic flow. The isothermal assumption eliminates the need for solving the 
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energy equation, but more complex assumptions can be used, e.g. a polytropic law. The system 
of equations to solve become mass conservation, momentum conservation, and the steady state 
induction equation (e.g. Pelletier & Pudritz (1992)). Respectively, these are 

V-(pu)=0 (2) 

pu • Vu = -J-(VxB)xB - VP - /)V$ (3) 
Vx(uxB) = (4) 

where <3? is the gravitational potential due to the central source. In cylindrically symmetric 
coordinates the physics of MCL disk winds can be cast in axisymmetric form where the velocity 
and magnetic fields are decomposed into toroidal (4>) and poloidal (r, z) components, 

B = B p + B^ (5) 
u = Up + = Up + Or <fi (6) 

where O is the rotational frequency of the plasma at a point (r, z). Axisymmetry and V • B = 
allow Bp to be expressed in terms of a magnetic flux function, a(r, z), such that 

Bp = -(Va x (7) 

The magnetic surfaces, generated by rotation of a poloidal field line about the axis, are surfaces of 
constant a(r, z). 

For axially symmetric configurations, (2) and (4), imply that the poloidal velocity is always 
parallel to magnetic surfaces, that is 

pup = k(a)B p , (8) 

where the function k is constant on a flux surface (KPOO equation 7). The induction equation and 
V • B = also imply 

pu^ = k(a)B (f> + prO G (a), (9) 

where G (a) is a constant that is approximately equal to the angular velocity in the disk where 
the magnetic surface is tied (KPOO equation 9). 

Using V • B = 0, Eqn. (2), and the azimuthal component of Eqn. (3), the conserved angular 
momentum per unit mass 1(a) becomes 

1(a) = Or 2 - = const(a) = £l r 2 A , (10) 

where the last equality follows from using (9) and finding the value at the "Alfven radius" r^(a) 
for a magnetic flux surface anchored to the disk at radius r (a) (e.g. KPOO equation 11). The 
Alfven radius r^(a) defines the radial coordinate of the point along a poloidal magnetic surface 
when the outflow speed on that surface equals the local poloidal Alfven velocity ua = B p / ^J-iirp. 
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The expression for conservation of angular momentum (10) thus relates contributions from the 
angular momentum of matter and magnetic torques to a fiducial value associated with the Alfven 
point. 

The poloidal component of Eqn. (3) can be integrated using (9) and (10) to give the total 
conserved specific energy U(a) carried by the wind on magnetic surfaces in both kinetic energy 
and Poynting flux. By defining w as the enthalpy per unit mass, the generalized Bernoulli integral 
then emerges (KP00 equation 14) as 

^(u 2 p + Q 2 r 2 ) + $ + w + Sl (Q r\ - Qr 2 ) = U{a) = const(a). (11) 



2.3. The Outflow Speed and Mass Outflow Rate 

We now assume a cold wind (such that w can be ignored). Because little acceleration occurs 
outside the Alfven surface we assume that 

Uoo = fUA (12) 

where / > 1. We will solve for ua and constrain / below. 

To calculate ua, we now solve the momentum equation for r < va in the rotating frame. For 
a cold wind, the dominant force components are the centrifugal force and the gravitational force. 
For r < TAi the magnetic field lines can be assumed to rotate rigidly with the angular speed of 
their foot points at r = ro- It is straightforward to see that the steady state radial momentum 
equation can then be written (Blandford & Payne (1982)) 

u .V^-^(^-^). (1 3) 

To keep the analysis as simple as possible, consider the initial launch to be highly inclined, such 
that z « r. Such an approximation is reasonable since (13) implies that the field inclination 
from the disk must make an angle < ir/4 to the disk plane for launch. This result for the field 
inclination follows from (13) by considering z « r and (r — ro) << ro and carrying out an 
expansion of the resulting force along the field line to second order in z 2 jr 2 and (r — ro ) 2 /r 2 . The 
reason for an outward force is the gain in centrifugal force at r > ro from the fact that "rigid" field 
lines enforce co-rotation. (Note in the following sections we will use the fact that the dominant 
contribution to the magnetic luminosity comes from the field lines anchored at the innermost 
radius r; L impling that ro = rj.) 

Integrating (13) along the field line for the case z « r and taking the result at r = ta gives 

4-^(1 + 24) =^Va, (14) 
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where we see that 

1< q 2 = (1 + 2rl/r\) < 3. (15) 

Now 

u\ = B 2 A /A7T PA , (16) 

where p A is the mass density at r A , and 



r\uA 47rr A fL q'' 



from the mass continuity equation, where £ is the solid angle < 4ir corresponding to the launched 
outflow mass loss rate M w j whereas M w = 4ttM W: i/^ is the effective mass loss rate were it 
quasi-spherical. 

Using (14), (16) and (17), we then obtain 

*- = H$T- <18> 



But using (12) and (1) we get a separate equation for M w , namely 



M - 2 B2aTa (V)\ 



By setting (18) equal to (19) we obtain f 2 q = 2. Thus from (15) we must have 

2 



3 1 / 2 



< r < 2. (20) 



This is a narrow range of / and for our crude order of magnitude estimates we will take 
/ = 1.2, a value right in the middle of the allowed range. For this choice of /, we obtain q = 1.4. 
For simplicity, we will use these values in what follows. From (15) these values imply 



r A = 1.27r . (21) 

Using these in (12) and (18) then give 

~ 2.in r . (22) 

and 

M w ~ 0.9^. (23) 
In the next subsection we obtain an expression for B A ■ 
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2.4. Magnetic Field Strength 

Magnetic fields may form in these disks via dynamo processes (Reyes- Ruiz & Stepinksi (1995), 
BFW01). The topology of such a field, (the ratio of poloidal B p and toroidal B§ in the disk) and 
its subsequent value in the coronae remains a subject of considerable discourse (e.g. Blackman 
2003). Here we assume that a field produced by dynamos can drive a disk wind in the manner 
described in the last section. This means that whatever combination of toroidal and poloidal 
field is produced in the disk, we assume a primarily polodial field in the corona where the wind 
launches. 

Eqn. (4) and V • B = imply that the radial magnetic field falls off with r 2 along the field 
line to ta- Thus 

B\ = B 2 (r /r A ) 4 = 0AB 2 , (24) 

where Bq is the poloidal field at the disk surface and we have used (21). The square of the surface 
poloidal field Bq can be estimated to be lower than the midplane poloidal field squared B\ by the 
density ratio to the 4/3 power (flux freezing). The density falls by a factor of e~ n , where n is the 
number of scale heights above the disk from where the wind is launched. We then have 

Bi = e-^Bj~4ne-^p d aid (25) 

where a ss is the disk viscosity parameter (Shakura & Sunyaev 1973), c s is the sound speed, and 
Pd is the disk density. The latter similarity in (25) follows from estimating the disk mean poloidal 
field from a helical mean field dynamo (Blackman, Frank & Welch (2001), Blackman 2003). From 
mass conservation in the disk, the disk density satisfies 

where ho is the disk scale height at r = ro and v r = a ss c s ho/rQ is the disk radial accretion velocity, 
and Md is the mass accretion rate. Combining this with (24), (25) and (26) then implies 

B 2 = ss = 2 5 £2 (27) 

n 



2.5. Compiling the Formulae for Application to Observations 

Combining (27) with (23) and using c s = Qho for an accretion disk, we obtain 

M w = 0.72e" 4n / 3 a ss ^M d (28) 

tlQ 

for the mass loss rate. Combining this with (22) and (1) gives 

L mag ~ 1.6e- 4n / 3 a ss ^M d n 2 rl (29) 
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The scaling for the momentum input II = M w Uoo of the wind can then also be estimated via 

n ~ L mag / Uoo ~ 0.76e- 4 ™/ 3 a ss ^M d n r . (30) 

ho 

In the above three relations, there is some freedom in choosing n, the number of scale heights 
above the midplane from which the wind launches (i.e. the location at which the plasma becomes 
magnetically dominated). One expects the corona to become magnetically dominated after one or 
two scale heights thus n = 1.5 is reasonable choice for the point at which we expect the ratio of 
thermal to magnetic pressure {0) to make the transition (3 > 1 to (3 < 1. Using this we obtain 

M w ~0.1a ss ^M d , (31) 

L rnag ~ 0.22a ss ^M d n 2 rl (32) 
no 

and 

fl~ L ma g/ Uoo ~0.Ua ss ^M d n r . (33) 

no 

The last two expressions give estimates of the rate that energy and momentum are input by 
MCL disk winds to the ambient medium, and (22) gives the asymptotic wind velocity 



Uqo ~ 2.m r . (34) 
These can be used for comparison with observations. 



3. Disk Accretion Rate in PNe 

In order to produce a more detailed comparison of MHD disk winds with PNe it is necessary 
to have a model for PNe accretion disks. In particular it is necessary to know the accretion rate 
M d , the inner disk radius as a function of time. 

It is unlikely that an accretion disk could survive the long main sequence lifetime of a PN 
central star. Thus, unlike YSOs and AGN, accretion disks in PNe systems must form via binary 
interactions Disks may form around secondaries via Roche lobe overflow or accretion of the 
dense AGB wind (Morris 1987, Mastrodemos & Morris 1998). Such systems would be similar to 
symbiotic stars (Corradi et a/.2000). Accretion disks could also form around the primary after 
CE evolution and disruption of the secondary star (Soker & Livio 1994, Soker 1998, Reyes-Ruiz & 
Lopez 1999). 
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Mastrodemos & Morris (1998) carried out detailed simulations of the first scenario. Using a 
SPH method they modelled the gravitational interaction of a dense AGB wind with a lower mass 
companion. They found steady accretion disks around a white dwarf companion orbiting a AGB 
star with M a gb ~ 10~ 5 M yr -1 . The ratio of Md/M a gb ~ -05 — .005 they found in their models 
is consistent with expectations from basic theory (Frank et al. (2002)) 

Accretion disks may also form via disruption of the secondary after CE evolution (Soker 
& Livio 1994). This model implies a finite lifetime for the disk as the mass reservoir of the 
disrupted companion is slowly drained onto the primary. A description of disk formation in PNe 
has been given in Reyes- Ruiz k, Lopez (1999). Envelope ejection occurs via transfer of angular 
momentum during which the secondary falls to a separation such that it catastrophically overflows 
its Roche lobe and forms a disc around the primary. There are a number of important constraints 
on the properties of binaries which would lead to disk formation in this way. Reyes-Ruiz &; 
Lopez (1999) found that systems with a primary consisting of an evolved AGB star with mass 
M* pa 2.6 — 3.6 M Q , a low mass secondary (< 0.08 M ) and an initial binary separation of 
< 200 i? may produce disks. The AGB star will shed most of its mass during the common 
envelope ejection, leaving a post-AGB stellar core surrounded by a thin shell. 

Reyes-Ruiz Sz Lopez (1999) find the disk accretion rate to evolve in time with in a power-law 
manner. 

• / t \ _5/4 

M d = M do [y^J M ©y r ( 35 ) 

(36) 

Typical values of the scale is Mdo = 10 -3 M Q yr -1 . 

4. Disk Winds Models for PNe and pPNe 

We wish to understand if disk wind models can account for outflows in PNe and pPNe and to 
characterize the parameters for the winds. The momentum II = M w u w and energy E = ^M w u\, 
injection rates for PNe are easily approximated. For "classic" PNe total mass of M pn pa .1 
Mq must be accelerated to velocities of u pn pa 40 km/s in a timescale of order At pn pa 10000 y 
(Note, here the subscript pn refers to observed properties of the total nebula which reflects material 
swept-up by the wind we model in this paper.) This gives II = M pn u pn / At pn pa 10 27 g cm/s 2 and 
E = M pn u 2 pn /At pn pa 10 34 erg/s. 

Recall that Bujarrabal et al. (2001)) found high total outflow momentum 
10 36 < U/(g cm s- 1 ) < 10 40 and total outflow energy 10 41 < E/(erg s' 1 ) < 10 47 in an 
extensive sample of pPNe. These values can be converted into momentum and energy injection 
rates using an assumed injection or "acceleration" timescale At, (II = U/At, L = Ej At. Recall 
that the values for II and E quoted above can't be explained via radiation wind driving. Note 
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also that there is uncertainty about At, the injection timescale, but most likely At < 10 3 y. The 
question which then arises is: can energy and momentum budgets be met with disk wind models? 

The results of section 2 demonstrate that disk wind solutions are sensitive to the location 
of footpoints for the flow r$ and the accretion rate M^. In what follows we use assume that the 
inner edge of the disk extends to the stellar surface and use rO = rj = r*. We shall see that the 
accretion rate is the parameter which becomes most important for obtaining solutions for pPNe 
since the requisite high outflow momenta will require high values of Md- One means of achieving 
high accretion rates will be to use the model described by RRL99 where accretion rates as high as 
a few times 10 -4 M yr _1 are possible for short periods as the disrupted companion's mass is fed 
onto the surface of the primary. 

PNe Solutions: We first consider the case of a "classic" PNe. In this case we would consider 
that the star which produces the jet is a proto-WD with an AGB companion (Soker & Rappaport 
2000). Thus accretion rates of ~ 10 -6 are resonable. To evaluate the expressions above we 
choose parameters for a canonical PN with mass M s = 0.6 M and a disk with a = 0.1 and 
ro/h = 10. 

If we assume typical PNe central star parameters (T* = 10 5 K, = 5000 L such that 
ri = 1.64 x 10 10 cm) we find the following conditions for the wind from equations 31 - 34 (note we 
use u w = Uoo), 

M, = 1 x I"' 7 Mo yr- ( 10 - 6 yr -l ) < 37 > 

- = L25xio3 ^»(^) i/2 fer /2 - < 3s > 

(39) 

Thus using typical conditions for PNe central stars, the scaling relations derived from the MHD 
equations yield disk wind parameters well matched with observations. 

pPNe Solutions: While the mass loss rates and velocities are known for PNe winds the 
situation for pPNe is not as clear. In general what is observed in pPNe is the total mass in the 
outflows. Mass loss rates must be inferred from the estimates of pPNe acceleration timescales. 
Velocities are also uncertain in the sense that the winds themselves may not be observed directly 
but only properties of swept-up material may be determined. 

We assume a model post- AGB star with mass M s = 0.6 M Q , a temperature of T* = 10, 000 K 
and a luminosity of = 5 x 10 3 L which, assuming a blackbody, yields a radius of 
r* = 1.6 x 10 12 cm = 23 R . Note that such a star has an escape velocity of u esc = 98 km s~ l 

From our previous discussion it is clear that achieving the high momentum input rates 
observed in pPNe via MCL disk wind models will necessitate high accretion rates. Thus we adopt 
an acceleration timescale of At = 200 yr and an accretion rate of Md = 1 x 10 -4 M yr -1 . This 
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value of Md is the 200 year average of that found by RRL99 for their case A (Mdo = 1.6 x 10 3 
M yr _1 , equation 35). Once again we choose a disk with a = 0.1 and r /h = 10. 

Assuming ro = r* along with the other parameter values given above, the key wind quantities 
M w , and L m = L w , E can again be determined from equations 31, 34 and 32, 

M, - 1. x 10- M yr- ( 10 _ 4 ^ yr _, ) (40) 

u w ~ 146 jfem/a ( ) I — ^- ) (41) 



.6M 7 V23R 



L m , a, X I0V. ^,^3,^ J (42) 

Note that L w At « 10 44 , a value in the middle of the range found by Burharrabal et al 2001. Note 
also that the solution above has Uqo 1.5u esc . Since « u esc the higher velocity outflows seen 
in some pPNe would require more disks around more compact central sources. We note here that 
the observed momenta and energy in the outflows comes primarily from swept up circumstellar 
material. Thus, as is the case in YSO molecular outflows (Bachiller 1996), the energy and 
momentum budget of the disk driven outflow must be sufficient to power the observed outflows 
via so-called prompt or shock driven entrainment. 

Given a model for the temporal history of the disk accretion, the total energy and momentum 
for the outflows can be found. Replacing in equation 35 with Md(t) from the relations derived 
by Reyes-Ruiz Sz Lopez (1999) and integrating with respect to t gives 

1/4 1 

(43) 

' i 

(44) 



E = J ^Mwulcdt w 1.3 x 10 44 erg 1 - f^Jll} ' 

n = /M.^i.sxio-^- 



i yr \ 

t ) 



These results show that the MCL disk wind models can achieve both energy and momentum 
injection rates as well as the total energy and momentum required to account for many pPNe 
described by Bujarrabal et al (2001). The total energy and momenta budgets we find from these 
solutions fall well within the range of pPNe outflows with momentum excesses. Taken together 
with our previous calculations for "classic" PNe winds these results confirm the predictions of 
Blackman, Frank & Welch (2001) that magnetized disk winds can account for much of the outflow 
phenomena associated collimated outflows in the late stage of stellar evolution. 

The results above indicate that collimated flows which form from transient disks in the pPNe 
stage will appear as dense knots in mature PNe flows. This may also serve to explain the presence 
of so-called FLIERS (Fast Low Ionization Emission Regions: Balick et aZ.1994) seen in some PNe. 
The mass loss rate in the winds derived above rapidly decrease with time. Thus the bulk of the 
jet's mass will lie near its head. As the material in the disk is accreted onto the star the jet will 
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eventually shut-off leaving the dense knot to continue its propagation through the surrounding 
slow wind. 

When the star makes its transition to a hot central star of a PNe its fast, tenuous spherical 
wind sweeps up a shell of the slow AGB wind material. The shell's expansion speed will typically 
be of order 40 km/ s and it will not catch up to the head of the jet. Thus during the PNe phase 
the jet head will appear as a dense, fast moving knot which should lie outside the PNe wind blown 
bubble. We note that masses of FLIERs are estimated to be of order 10~ 4 — 1CT 5 M which is 
reasonable for the models presented above. FLIER velocities can be lower than the « 100 km/s 
calculated above but deceleration of the jet head will occur via interaction with the environment. 
We note also that hydrodynamic simulations of PNe jets in which the jet ram pressure decreases in 
time (as would occur for our model) show characteristic patterns of backward pointing bow-shocks 
(apex pointing back towards to the star, Steffan & Lopez 1998). If such results are robust, the jets 
produced by disk winds in our scenario above may also yield similar morphologies. 



5. Discussion and Conclusions 

While purely hydrodynamic models for PNe shaping provide an adequate description for many 
large scale features seen in purely morphological studies, jets and hourglass shaped bipolar nebulae 
appear to strain their explanatory power. In addition, the momentum and energy associated with 
many pPNe appear to be orders of magnitude larger than can be accounted for with radiation 
driving from the central star (even when multiple scattering is taken into account). In this paper 
we have derived and applied scaling relations derived from time-independent axisymmetric MHD 
equations to the winds and wind-driven outflows in PNe and pPNe. 

For "classic" PNe (a central star with % » 3 x 10 4 K) we find magneto-centrifugal winds 
can account for typical observed wind properties. We find only a 1% efficiency of accretion of an 
AGB companion wind is required to produce reasonable results. 

Our results for pPNe show that momentum excesses need not occur for outflows driven 
by MCL winds. While this is encouraging in terms of finding a mechanism for driving pPNe 
outflows the solutions require fairly high accretion rates (> 10~ 5 M Q yr" 1 ). It is not clear if 
such conditions can be achieved with the frequency required by observations. While solutions of 
Reyes-Ruiz &; Lopez (1999) yield accretion rates and time dependencies which lead to the correct 
outflow momenta and energetics, their models place fairly stringent limitations on the nature of 
the binaries that form disks from disrupted companions. If accretion onto undetected compact 
orbiting companions is invoked (Soker & Rappaport (2001)) then higher values of Md may not be 
required. 

It is worth noting that these models imply an outflow from both the AGB wind and the jet. 
The AGB wind can, in princple be sculpted by its own magnetic forces or can be swept at later 
times when a radiation-driven wind from the exposed core is initiated. Thus this class models 
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implies the possible mis- alignment between the jet and the main body of the nebula. In Blackman, 
Frank & Welch (2001) a model for such multi-polar outflows was proposed in which magnetically 
launched outflows are driven from the both the disk and a non-aligned rapidly rotating AGB core. 
In this paper we did not consider the AGB wind to be magnetized and considered the field in the 
disk to be arise via a turbulent dynamo there. If the AGB wind is magnetized then shaping could 
occur with an alignment that is uncorrelated with that of the collimated disk wind. We note that 
current MHD models of AGB winds however tend to support the development of a torus (Matt 
et al 2001) rather than a jet. Jet formation via initially weak fields in post- AGB winds have been 
explored (Garcia-Segura et al 1999) however these models require higher winds speeds than occur 
on the AGB. Finally it is worth noting that recent simulation results by Matt, Blackman & Frank 
2004 confirm that the exposed rapdily rotating magnetic core model can produce well collimated 
outflows. Regardless of these cavets the point remains that collimated disk winds and AGB wind 
systems can, in principle, explain multi-polar outflows. 

We note that the presence of molecules in fast moving outflows driven by winds, as proposed 
in our studies, is not new. Studies of YSO molecular outflows have revealed "fast molecular" 
material moving at speeds of more than 25 km/s, which is the nominal dissociation speed for H2 in 
J-shocks. A number of theoretical proposals to resolve the issue currently exist including magnetic 
precursors in J-shocks (Harigan et al 1989), high velocity C-shocks (Smith & Brand 1990), 
accelerating shocks (Lim 2003) and the presence of clumped or inhomogeneous gas (Hartquist & 
Dyson 1987). Thus the behavior we see in proto-Planetary nebulae also exists in YSOs where 
compelling evidence for MHD driven disk winds already exists. The question of how molecules 
can survive after being swept-up in high velocity flows remains unanswered for pPNe and YSOs 
but is currently an active area of research with a number of competing models currently under 
investigation. 

We note that a robust prediction of our models is the ratio of wind mass loss rate to accretion 
rate, i.e. M w /M a f» .1. This is true for most MCL disk wind models and can be seen as a target 
prediction which can be explored observationally. 

Thus the question which must now be addressed is can conditions for either high accretion 
rate disks in pPNe or symbiotic type accreting companions be made to embrace enough systems 
to account for the statistics of Bujarrabal et a/.(2001). There are many uncertainties concerning 
the formation of disks via the disruption of the companion after a CE phase. The disruption of a 
secondary may not really lead to an "accretion disk" around the primary due to both tidal and 
hydrodynamical disruption. The requirement that material from the disrupted companion "push" 
trough the envelope may not lead to a disk but rather a dense and clumpy expanding torus. 
A detailed discussion of the process however remains outside the scope of our paper. Readers 
wishing to consider the viablity of these models are encouraged to review Soker & Livio (1994) 
and Ruiz-Reyes & Lopez (1999). 

There is an important difference between the two disk formation scenarios discussed here. In 
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the case where AGB wind material is captured to form a disk, one expects the abundances of the 
jet and AGB outflows to be the same. In the disrupted companion scenario the disk/jet will have 
formed from a different star in a different evolutionary state. Thus we would expect abundance 
differences between the AGB and jet outflows. It may, therefore, be worthwhile for observers to 
search for abundance gradients between different components of multipolar flows. 

We note finally that this paper comprises a step beyond Blackman, Frank & Welch (2001) 
in establishing the efficacy of MHD paradigms for pPNe/PNe in which strong magnetic fields 
play a role in both launching and collimating the flows. Future models should aempt to include 
more detailed description of the physics of MCL disk launching in these systems including time 
dependent models, (see for example (von Rekowski 2004) 
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